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Abstract 



We formulate a consistent model predicting, in the leading order approximation, 
r"! I maximal atmospheric mixing angle, vanishing reactor angle and tan0i2 = 1/0 where 

(j) = {1 + \/5)/2 is the Golden Ratio. The model is based on the flavour symmetry 
' ^5 X Z5 X Z3, spontaneously broken by a set of flavon fields. By minimizing the scalar 

• potential of the theory up to the next-to-leading order in the symmetry breaking 

On . parameter, we demonstrate that this mixing pattern is naturally achieved in a finite 



portion of the parameter space, through the vacuum alignment of the flavon fields. 
The leading order approximation is stable against higher-order corrections. We also 
compare our construction to other models based on discrete symmetry groups. 
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1 Introduction 



Neutrino oscillations have revealed unexpected features of the flavour problem. The mixing 
matrix in the lepton sector is completely different from that in the quark sector and the 
relative hierarchy between neutrino masses, given in terms of the ratio between the two 
independent squared mass differences, is small compared to the hierarchy for charged 
fermions. The regular pattern observed in the quark sector strongly suggests that quark 
masses and mixing angles should be explained by some dynamical principle which, in a 
unified picture, should eventually apply also to leptons. In searching for such a common 
principle a valuable tool is that of spontaneously broken flavour symmetries, which might 
be tentatively identified from the present data. 

The lepton mixing matrix Upmns still suffers from large uncertainties. The parameters 
related to CP violations are totally unknown at present. The reactor angle ^13 is the 
smallest mixing angle, but there is only an upper bound on it and its value can range from 
zero to about 0.2. The atmospheric mixing angle 623 is compatible with being maximal, but 
deviations from maximality are still allowed to some extent. The most precisely measured 
angle is the solar angle 612, which is large but not maximal, with a Icr uncertainty of less 
than 2 degrees: 

sin2^i2 = 0.304iH?^ [I] , sin2^i2 = 0.32ll0:0i [2] . (1) 

Despite the lack of a precise knowledge of Upmns, the present data can guide us in 
searching for a first-order approximation providing the basis of a theoretical description. 
Spontaneously broken flavour symmetries are suitable to this type of description since 
they typically generate a mixing matrix at the leading order (LO) of the expansion in the 
symmetry breaking parameter, with small sub-leading corrections originating from the 
higher-orders. In the quark sector a natural first approximation is given by Vckm = 1- In 
the lepton sector the identification of such a first-order approximation is less obvious and 
several candidates have been proposed so-far. Most of them have in common atmospheric 
and reactor angles, fixed to the values sin^^^23 = 1/2 and 6*13 = 0, and they only differ in 
the solar mixing angle. Tri-bimaximal mixing (TB) is perhaps the most studied pattern 
[5]. It predicts sin^^i2 = 1/3, which is within two standard deviations from the current 
best value. TB mixing can be reproduced at the LO in many models based on discrete and 
continuous flavour symmetries [1]. A minimal construction is based on A^. In bimaximal 
mixing (BM) the solar angle is maximal, sin^6'i2 = 1/2, outside the presently allowed 
range O |6]. To reconcile the LO approximation with the data, the expansion parameter 
should be not-too-small, of the order of the Cabibbo angle. Sub-leading corrections of this 
size are expected to affect also other parameters, such as 613, which is thus predicted close 
to the present experimental upper bound. 

Another plausible mixing pattern is the one where sin^ ^23 = 1/2, ^13 = and tan 612 = 
where (j) = {1 + V5)/2 is the Golden Ratio (GR) This pattern, called GR the 

hereafter, is the focus of the present paper. We have 

sin^ 9i2 = = — ^ ^ 0.276 . (2) 

5 + ^5 
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This value is about two standard deviations below the experimental range and can be 
brought inside the allowed interval by a small NLO correction, of order 0.05 radiants. 
The GR characterizes several properties of the icosahedron and a natural candidate for 
the flavour symmetry giving rise to (j2]) is the icosahedral one, related to the group A^. 
This relation was pointed out for the first time while trying to connect the value of the 
solar angle to the Golden Ratio [8]. Indeed there have been attempts to construct a model 
based on the symmetry [11] to the purpose of reproducing the GR mixing pattern, 
but a complete model does not exist to date. Recently, the group A^ was also applied 
to a scenario with a fourth lepton family [12], while the double cover of the icosahedral 
group was used to reproduce the quark mixing [13]. In models based on spontaneously 
broken flavour symmetries a crucial feature is the discussion of the vacuum alignment. The 
family symmetry is broken by the VEVs of fiavon fields and the desired mixing pattern is 
intimately related to the directions of these VEVs in flavour space. In a complete model 
the VEV alignment should occur naturally, as the outcome of the minimization of the 
energy density of the theory. To our knowledge none of the existing proposals of what we 
called the GR pattern have solved the vacuum alignment problem. 

Aim of our work is to build a complete model based on the family group A^ and 
reproducing the GR mixing pattern in a natural way. In section 2 we will discuss how the 
symmetry group A^ can be used to generate the GR mixing pattern. We show that the 
most general neutrino mass matrix giving rise to the GR mixing pattern is invariant not 
only under the fi — t exchange symmetry, but also under a parity transformation 5*. At 
the same time the requirement of a diagonal charged lepton mass matrix is guaranteed 
by the invariance under a transformation T of order five. We will see that the elements 
5* and T actually generate the group A^. Some properties of this group will be recalled 
in section 3. In section 4 we define our model by assigning all fields to representations of 
the flavour group. We also list the fiavon VEVs needed in order to accomplish the desired 
symmetry breaking and generate the GR ratio. The minimization of the scalar potential 
is performed in section 5 and 6, respectively at the LO and NLO. In section 6 we also 
enumerate the other sources of corrections to the LO approximation. We compare our 
model to other proposals in section 7 and then we conclude. 



2 A family symmetry for the Golden Ratio 

We start by analyzing the property of the most general neutrino mass matrix leading to 
the Golden Ratio (GR) prediction for the solar mixing angle. We chose a basis where the 
mass matrix for the charged leptons m/ is diagonal. More precisely, it is sufficient that the 
combination m\mi is diagonal, so that there is no contribution to the lepton mixing from 
the charged lepton sector. We should also make a choice for ^23 and ^13. To begin with we 
assume a leading order approximation where sin^^23 = 1/2, sin^6'i3 = and tan 6*12 = 1/0 
where = (1 + \/h)/2 is the GR. We look for the most general neutrino mass matrix 

^An alternative proposal f5I relating the Golden Ratio to the lepton mixing assumes cos0i2 — (/)/2. 
Consequently we have sin'^ 612 = -1(3 — </>) ~ 0.345, about two standard deviations above the experimental 
value. In |10) this prediction was deduced from the symmetry of the dihedral group -Dio- 
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rrii, leading to this mixing pattern. Such a matrix can be constructed by acting with the 
corresponding mixing matrix Ugr on a generic diagonal neutrino mass matrix: 



Ugr diag(mi,m2,m3) Ul 



GR 



(3) 



In a particular phase convention, the matrix Ugr representing our mixing pattern is given 
by: 



GR 





COS 6*12 


sin 6^12 


\ 




sin 612 


cos 012 


1 




V2 


V2 






sin 6^12 


cos 6*12 


1 


\ 


V2 


V2 





(4) 



with tan^ 



12 



). By applying eq. ([3]) we find a matrix of the form: 




(5) 



with coefficients x, y, z and w satisfying the following relation: 

z + w = X — \/2y . 



(6) 



The matrix in eq. ([5]) is the most general one giving rise to ^13 = and 623 maximal. The 
constraint of eq. (E]) arises from further specifying the solar mixing angle. 

The matrix in eqs. (IMSl) can be completely characterized by a simple symmetry re- 
quirement. Indeed, it is invariant under the action of the two unitary transformations: 



U 




s 




(7) 



which satisfy 



[S, U] = 



and generate a group Gj, = Z2 x Z2. Conversely, the requirement of invariance under U 
and S completely characterize rriiy in eqs. (ISM]). Namely, given a generic neutrino mass 
matrix rriiy the most general solution to the equations: 



S'^ nil, S 



m,. 



(9) 



with U and S given in eq. (^^, is the mass matrix defined by eqs. ([5]) and (EI). 

In the chosen basis, where m\mi is diagonal, there is no contribution to the lepton 
mixing from the charged lepton sector and the mixing matrix Ugr originates only from 
the diagonalization of rriy. To construct a model for the desired mixing pattern, we should 
require that a diagonal m\mi arises naturally, as the general solution of a symmetry or 
dynamical requirement. For instance, we can require that the charged lepton sector is 
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invariant under a family group Gi enforcing a diagonal mjmi. In our LO approximation 
the groups G^, and Gi should be seen as the residual vacuum symmetries characterizing the 
neutrino sector and the charged lepton sector, respectively. Such a configuration can be 
induced by the spontaneous breaking of some family symmetry Gf, through the vacuum 
expectation values of two different sets of flavons that selectively couple to neutrinos 
and to charged leptons. It is not strictly necessary that Gf entirely contains Gi and Gi^ 
as subgroups, since a part of the residual symmetries can arise accidentally, due to the 
specific field content of the model, as the baryon and the lepton numbers arise as accidental 
classical symmetries in the standard model. A natural candidate for the family symmetry 
Gf giving rise to the GR prediction for the solar mixing angle is the proper symmetry 
group of the icosahedral, the alternating group [2]. One of the possible presentations 
of is in term of two generators S and T satisfying: 

= {STf = 1 and = 1 . (10) 

We make the following ansatz: we identify the matrix 5* in eq. ([7]) with the generator S 
of A^. Given the explicit form of the generator S, the algebraic relation ffTOj) allows to 
determine the matrix corresponding to the generator T. We find: 



/ 1 








\ 




27ri 









e~ 











Sni 




Vo 







/ 



This is an encouraging result. Indeed the condition 

{ra\mi) T = (m|mi) (12) 

requires m\mi to be a diagonal matrix and the natural candidate for the subgroup Gi is 
the group generated by T. We look for a model invariant under the family symmetry 
A5, where, after spontaneous breaking, the residual symmetries of the neutrino sector and 
of the charged lepton sector are those generated by {S, U) and T, respectively. It turns 
out that, up to an irrelevant overall sign, the /i — r symmetry is an element of A5, that can 
be expressed in terms of the generators S and T of eqs. ([7]) and (ITT]) §. By construction 
the model predicts GR for the solar mixing angle. We expect that higher order corrections 
produce small deviations from the LO predictions = and 6*23 = tc/A and tan^^i2 = 1/0. 

Notice that this approach automatically guarantees the independence of the mixing 
matrix Ugr and the other physical results from the base choice. Indeed, in a generic 
basis where the generators are Xq = Q X Q'' {X = S,T,U), Q denoting a unitary 3 x 
3 matrix, in general the combination mlmi is no more diagonal and the neutrino mass 
matrix m,^ have a texture different from the one in eqs. ([MSD. However, as a result of the 
residual symmetries, m|m/ is diagonalized by Q, whereas rriu is diagonalized by {Q Ugr), 
the physical mixing matrix remaining Ugr- 

^We thank Claudia Hagedorn for pointing this out to us. 
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3 The group 



The group is the group of the even permutations of five objects. It is the proper 
symmetry group of two of the five Platonic sohds, the icosahedron and the dodecahedron. 
It has 60 elements that can be grouped into five conjugacy classes with 1, 12, 12, 15 and 20 
elements. The five irreducible representations are the invariant singlet, two inequivalent 
triplets, a tetraplet and a pentaplet. The characters of A5 are collected in Table 1. The 



A. 






12C^l 


12Cf 


15Cf 


20Cfl 

5 




1 


1 


1 


1 


1 




3 




(1-0) 


-1 







3 


(1-0) 





-1 







4 


-1 


-1 





1 




5 








1 


-1 



Table 1: Characters of the A^ group. 



products of two A^ representations can be decomposed according to the following rules 



3C 


?)3 = 


(1 + 5)5 + 3^ 


3'g 


)3' = 


(1 + 5)5 + 3:4 


3(2 


)3' = 


4 + 5 


3C 


54 = 


3' + 4 + 5 


3'C 


54 = 


3 + 4 + 5 


3(S 


5 5 = 


3 + 3' + 4 + 5 


3'C 


55 = 


3 + 3' + 4 + 5 


4C 


54 = 


(1 + 4 + 5)5 + (3 + 30a 


4C 


55 = 


3 + 3' + 4 + 5 + 5 


5C 


5 5 = 


(1 + 4 + 5 + 5)5 + (3 + 3' + 4)a 



where the suffices S{A) denote the symmetric(antisymmetric) property of the correspond- 
ing representation. The product between the singlet and any representation r gives r. 
As recalled in the previous section, A^ is generated by two elements S and T, with the 
presentation 

S'^^{STf^l and 7^ = 1 . (14) 
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The element S belongs to the class C[ ^ and the element T to the class C2 . We find useful 
to work in a basis where the T generator for the various representations is always diagonal. 
Since this choice is unconventional we list in Table 2 the matrices associated to S and T 
in our basis. In Appendix B we made connection with other basis used in the literature. 





S 




3 




' 1 V2 x/2\ 




diag(0, 1, 4) 


3' 


1 


/ 1 x/2 


/ 


diag(0,2,3) 


4 


f (0-3) (0 + 2) -v^ \ 

1 (0-3) (0 + 2) 
5 (0 + 2) V5 (0-3) 
\ -V5 (0 + 2) (0-3) -V^ I 


diag(l,2,3,4) 


5 


/ -1 ye -\/6 -^6 \ 

2-0 20 2(1-0) -1-0 
\ -v^ 20 1 + 2-0 2(-l + 0) 
2(1-0) 2-0 1 + -20 
-1-0 -2(1 -0) -20 2-0 ) 


diag(0,l,2,3,4) 



Table 2: S and T generators of in the basis where T is diagonal. 



Notice that the matrices S and T of the previous section coincide with those of the 
generators in the representation 3. We have derived the Clebsh-Gordan coefficients 

entering the decomposition of the representation products. They are given in Appendix 

A. 

4 A model with family symmetry 

In this section we define our model. We focus on the lepton sector and, to facihtate the 
task related to the vacuum alignment, we consider a supersymmetric model in the limit 
of exact supersymmetry (SUSY). SUSY breaking effects do not affect lepton masses and 
mixing angles. Among the fields in the lepton sector we include three gauge singlets and 
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the neutrino masses will be dominated by the contribution of a type I see-saw mechanism 
[T5] . A version of the model without see-saw, where the neutrino masses are described 
by effective higher-dimensional operators, is equally possible. It would lead to the same 
predictions for the lepton mixing angles. 

To start with we assign both the SU(2) lepton doublets / and the right-handed neutrinos 
v'' to the representation 3 of A5. We take the SU(2) singlets e*^, and as invariant 

singlets. Higgs doublets ilu,d are also singlets of ^5. In the neutrino sector we can 
write a renormalizable Yukawa coupling of the type [yH^H^^ the notation (...) standing 
for the combination of the fields in parenthesis giving an singlet. The product v'^v^ is 
symmetric and contains a singlet and a pentaplet of and, to discuss the most general 
case, we introduce two flavon chiral multiplets ^ and ^pg transforming as 1 and 5 of A5, 
respectively. They are completely neutral under the gauge interactions. In the charged 
lepton sector renormalizable Yukawa interactions are not allowed and we need additional 
flavons transforming as 3 under ^5. To solve the vacuum alignment problem a minimum 
of three triplets is needed, since trilinear interaction terms depending on less than three 
triplets vanish by the A^ symmetry. We include three triplets <y9, ^p' and neutral 
under the gauge interactions. An additional flavon singlet of ^5, is also introduced to 
implement the desired vacuum alignment. To avoid couplings of the flavon multiplets to 
the wrong sector we also need to enlarge the flavour symmetry. This is done by considering 
the group = ^5 x x Z3. In Table 3 we collect the chiral supermultiplets and their 
transformation properties under Gf. Notice that, at variance with other constructions 
based on flavour symmetries, we do not introduce the so-called driving fields. 





6^= 


If 




/ 




Hu,d 














As 


1 


1 


1 


3 


3 


1 


3 


3 


3 


5 


1 


1 


^5 





4 


1 














4 


1 











^3 


1 


1 


1 


2 


1 














1 


1 


2 



Table 3: Chiral multiplets and their transformation properties. 



The additional symmetry Z3 is a discrete version of the total lepton number and is 
broken by the VEVs of the flavons of the neutrino sector, ips, ^ and This symmetry 
prevents a direct mass term for u'^. The presence of the new factor forces each of the 
lepton multiplets e'^, and to couple to only one of the triplets ip' and y?", at the LO. 
The additional factors and Z^ play also an important role both in the construction of 
the flavon scalar potential and in the classification of NLO corrections. The superpotential 
for the lepton multiplets reads: 



w = y{i^'l)Hu+ \l-yi^iiy'iy') + y5i^si^'iy') 

+ yee'(^l)H, + y^i,-(^l)Ha + yrT''(^l)H, + ... 



(15) 



where dots stand for higher order operators and A denotes the cut-off scale. Notice that 
the LO Yukawa couphngs of the charged fermions are described by non-renormahzable 
operators. As we will see in section 5, where we will discuss the vacuum alignment, at the 
LO the flavons v?5, ^, V', ^' and acquire VEVs of the type: 





(-Y^(p + ?),-p,9,g,p) A 


(0 = 


s A 


(y^) = 


(m,0,0) a 


(^') = 


(0,m',0) a 




(0, 0, u") A . 



(16) 

where p, q, s, u, u', u" are dimensionless coefficients. Such a pattern completely specifies 
lepton masses and mixing angles, at the LO. Plugging the VEVs of ips, ^, ^p, ^' and y?" 
into the superpotential w and working out the invariant combinations, with the help 
of the results of the previous section and those of the Appendix A, we can find the LO 
mass matrices mi and vtiy. 

In the charged lepton sector, after breaking of A5, the relevant part of the superpoten- 
tial becomes 

yeU e^lelid + VpU' [l^l^xHd + VrU T^'lrHd . (17) 

There is no contribution to the lepton mixing from this sector and charged lepton masses 
are 

me = VeUVd , = y^u'vd , = Vru'vd , (18) 

Vd being the VEV of the neutral component of Hd- We might be surprised by the fact that 
mi is diagonal, since only the VEV of ip leaves the Z"^ subgroup generated by T invariant, 
while the VEVs of y?' and if" break Zj. We can understand this result by recalling that 
the flavour symmetry Gf contains a factor Z5, distinct from Z'^ . The VEVs of (p, ip' and 
ip" break x Z5 down to the diagonal subgroup Z^ contained in the product Zj x Z5. 
It is this residual group that guarantees a diagonal mi in our construction. 

Similarly, in the neutrino sector we read from eq. ( ITSjl the mass matrices M for the 
right-handed neutrinos and m^, the Dirac one: 



M = \/6 



yis 



1/5 



\yb{p + q) 
p_ 

p_ 

V2 



rriD 



2/5 



P 



V2 

yis- lybip + q) 



Z/5 



p 

V2 



\ 



yis- -y5{p + q) 
ym 



A 



(19) 




(20) 



where is the VEV of the Higgs doublet if^. Since M is fi — t symmetric and m^ is 
proportional to the matrix U in eq. ([7]), from the see-saw formula we have 



m,. 



m^M~ 



'mn = y^vlM 



(21) 
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We notice that M has precisely the structure given in eqs. (E]) and ^ and therefore both 
M and its inverse are diagonahzed by the mixing matrix in eq. ([1]) with tan = 1/0: 



Ugr Ugr = diag(mi,m2,m3) . (22) 

Therefore Ugr represents the contribution to the lepton mixing coming from the neu- 
trino sector, as desired. This resuh crucially depends on the VEV of the flavon penta- 
plet (fis, which will be derived from the minimization of the scalar potential in section 
5. We observe that such a VEV is left invariant by the action of the generator S, as 
can be immediately checked by multiplying the 5x5 matrix S of Table 2 and the vector 
-P, Q, QiP)- This is the reason why the residual symmetry of the neutrino 
sector contains the parity subgroup generated by 5". The presence of the /x — r symmetry 
is more subtle. Indeed the generator S of the 5 representation has three eigenvalues equal 
to one and the corresponding eigenvector can be parametrized as 

(-a/2/3(p + g + rcf)), -p + r, g + 2r(j), q,p + r) . (23) 

This is the most general VEV of ips that leaves S unbroken. It is easy to construct the 
corresponding neutrino mass matrix rriu and check that in the general case, with r 7^ 0, 
nil, is not yU — r symmetric. In our model it is the minimization of the scalar potential that 
selects the vacuum with r = 0, thus enforcing the yU — t symmetry. We will demonstrate 
this result in section 5. 

The charged fermion masses depend on three sets of independent parameters, which 
do not display a manifest relative hierarchy. It is easy to induce the correct hierarchy by 
assigning Froggatt-Nielsen U(l)i? charges 2q and q to e'^, and fj,'^, respectively [I6]. The 
spontaneous breaking of such U(l) p by the VEV of a scalar fields carrying a negative units 
of F explains why ye « « y^. In the LO approximation the spectrum of the light 
neutrinos is 

'^' = A + l + C ^ '^'= A + l-C ^ '^' = ^aTw- 
where A, B and C are complex parameters defined as 

^ = 6^1^^^^' B = y,{p + Aq)^J^, C = 3V5y,p-J^. (25) 

There are no special relations between the three complex parameters and thus we have no 
prediction on the neutrino spectrum, that can have both normal and inverted ordering. A 
moderate tuning among the parameters is needed in order to reproduce the ratio between 
solar and atmospheric squared mass differences. The mass combination entering neutrino- 
less double-beta decay, rriee, is given by |/^w^(M~^)ii and depends on the same parameters 
A, B and C. By expressing the absolute values of A, B and C in terms of Am^oi, Am^^^, 
the smallest neutrino mass and the phases of A, B and C, we can derive a range for |mee| 
as a function of the smallest neutrino mass by varying the available phases. We plot the 
result in figure 1, where the lighter region refer to the unconstrained case [17], and the 
darker one corresponds to the present model, in the LO approximation. 
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Figure 1: Predictions for \mee\ versus the lightest neutrino mass, for normal (red) and 
inverted (blue) mass ordering. In lighter regions the mass parameters are completely free. 
In darker regions they are correlated as prescribed by the LO approximation of our model. 

The dominant contribution from local effective operators to the light neutrino masses 
is (UhuhuC)- This operator is suppressed compared to the see-saw contribution, since the 
former is of order VEV/A"^, while the latter is of order 1/VEV. Moreover, it is easy to see 
that the new effective term preserves the /i — r symmetry and the Golden Ratio prediction. 

5 Vacuum Alignment 

The results of the previous section crucially depend on the assumed set of VEVs, eq. ( IT6|) . 
and the purpose of this section is to show that they derive from the minimization of the 
scalar potential of the theory, without ad-hoc tuning of the parameters involved. The 
transformation properties of the flavon fields allow to write the following contribution to 
the superpotential 



where dots stand for higher dimensional terms, which will contribute at the NLO. There 
are two independent cubic invariants that can be built from a pentaplet and they are 
denoted by the suffices 1 and 2 in Wf. There are no driving fields in our construction and 
the minima are derived by analyzing the F-terms of the flavons themselves. At the LO 
there is no mixing between ^, ips, that control the neutrino mass terms and (p, (f' and 
(f", that give rise to the charged lepton Yukawas. We can separately discuss the two sets 




(26) 
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of minima. We start from the neutrino sector. The condition 



is solved by 



-Mo^ + g,e = (27) 



Another set of conditions is given by 

^"^-0 . (29) 



To solve these equations it is convenient to move to the so-called Cummins-Patera basis 
flEi [T9] for the generators S and T. In this basis the generator S for the five-dimensional 
representation is diagonal, Sep = diag(+l, — 1, — 1, -|-1, -|-1). The exphcit form of T for 
the 5 representation in the Cummins-Patera basis as well as the unitary matrix relating the 
two basis is given in Appendix B. We denote the components of in the Cummins-Patera 
basis by 

ifs = iX,,X2,X,,Z,Z) , (30) 

where Xi {i = 1,2,3), Z and Z should be seen as independent complex quantities. The 
terms of the superpotential wj that depend on (ps are explicitly given by: 

Wf = g^i{Xl + Xl + Xl + 2ZZ) 

+ g2 [Z^ -Z^ - ^{Xl + u'^Xl + ujXI)Z + ^{Xl + uXl + w^X|)Z) 

+ ^3 {Z^ + Z^ + {Xl + u'^Xl + uXl)Z + [Xl + uXl + u'^XDZ 

- 4X1X2X3) + ... (31) 

where 

27rz 

UJ = e~ . (32) 

We have explicitly solved the equations (!29|l in this basis. We found no non-trivial solutions 
invariant under T and seven independent solutions invariant under S. They can be grouped 
in two pairs and a triplet. Each of these sets is closed under the action of the generator 
T. One of the triplet of solutions is given by: 

Xi=X2 = X3 = , (33) 

2gig^ 1 12 

~3Mo {g2-gzy'\92 + 9z?/'^ ' 

3Mo {g2-9,Y''{92 + 9,Y/'^ ' ^ ^ 

where we have also made use of eq. (128|) . The other two solutions belonging to the triplet 
are obtained by multiplying Z by uj{uj'^) and Z by uo'^iu). The condition X2 = X3 = 
correspond to the invariance under whereas Xi = is an additional specific feature of 
this set of solutions. 
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In each of these minima we have the basis independent result: 

89bl 1 



9M2 (g^-g^)(g^ + g^ 



(35) 



We have a finite multiphcity of minima and we chose the minimum in eqs. fl33|34p . Before 
coming back to our basis, we analyze the equation 



dwf 



(36) 



which, by making use of 

Mo 



and becomes 

9(92- g3)ig2 + gs) 



g^gl 



j/4 



Ml 







(37) 



This equation has non-vanishing solutions for which make non-trivial the solutions fl28|l 
and (l33l With the help of the unitary transformation relating the Cummins-Patera 
basis to ours we find that in our basis the minimum (l33l EH) translates into 



(3^ 



with 



1 



P 



2V2A 
1 

2y2A 



— I 



+ t z 



z 



(39) 



with Z and Z given by eq. (IMIl . We have recovered the pattern displayed in section 4. 
Notice that this result does not depend on the specific value of the Z and Z components, 
but rather on the conditions Xi = 0. In particular, the n — t parity symmetry is related 
to the vanishing of the Xi component. 

Moving to the fiavons ip, ip' and ip", the relevant part of the superpotential is given by 



Wf = Mi{(p'^) + M2{(p'(p") + g{(p(p'(p") + ... 

= Ml {(fl + 2(p2(p3) + M2 {(p[(Pi + ip'^vl + ) 

+ gi^i^'^'^l + + ^z^'x^l - ^i^W2 - ^2vWi - V'ivWi) + ••• 

The minima in the y?, y?' and (p" can be found by solving the system of equations: 

dWf _ Q ^"^Z _ n 







which, using a vectorial notation, can be written as 

2Mi + g (f' Xip" = , 2M2 ip" -g ipxip" = 







(40) 



(41) 



2M2 ip' + g ipxip' = , (42) 
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where x denotes the external product and, for any vector v = (fi,t'2,f3), we set v = 
{vi,V3,V2). To solve this system it is useful to recognize that the LO part of the su- 
perpotential that depends only on the fields if, if' and ip" is invariant under the linear 
trasformation 



(f — )■ QR Q V ; 

where y is a complex dimensionless parameter, 

/l 



1 



1 '> 



y 



(43) 



n 



\ 



1 -i 

- - , 

V V2 V2 J 



(44) 



and i? is a general complex orthogonal matrix, R^R = 1, depending on three complex 
parameters. As we shall see this invariance is accidental and is broken by the NLO 
contributions to the superpotential. By exploiting such an invariance we can always reach 
the particular minimum with = (ps = 0. It is easy to see that one such solution is given 

bjil 



^0 



If 



-^^ (1.0.0) 

9 



V2MIM2 

9 

V2M1M2 
9 



(0,1,0) 



(0,0,1) 



The general solution of the system (1421) is given by 



(p = flR Q ^{pQ 



1 

y 



QR Q 



-1 / 

^0 



(45) 



(46) 



The degeneracies related to R and y are accidental. Indeed the transformations of eq. 
fll5]) are not symmetries of our system, but rather accidental symmetries of the LO ap- 
proximation. It is easy to see that the symmetry related to the rotation R is removed al 
the NLO, by including operators of dimension four in Wf. As we will discuss in the next 
section, the inclusion of the most general set of operators of dimension four in the flavon 
fields, leads to the result R = 1, thus justifying the choice of the previous section. The 
symmetry under the rescaling y is removed by adding invariant operators of dimension 
five. 



6 Higher-order corrections 
6.1 Vacuum alignment 

The vacuum alignment discussed in the previous section is modified by the contribution to 
the superpotential from higher dimensional operators. If we denote by VEV the typical 

•^There is also another solution where the entries of (p' and ip" are exchanged. 
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vacuum expectation value of the flavon fields, we expect corrections to the LO minima of 
order VEV/K. These corrections can be kept small by asking VEV/ ^ 1. Nevertheless 
they play an important role in removing some of the degeneracy that affect the LO result. 
In the following discussion we include all NLO operators, that is operators of dimension 
four depending on the flavon fields. A complete set of invariants under the flavour group 
is given by 



^ 1 


( LOLO^ ( LOLO^ 


On 






vrH^ ) VrH^ ) 
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Q3 = 


ivVWv") 


Ql3 


= {{{^3^3)5^^3)3^) 


Qa = 


(y,V')(^'V") 


QlA 


= {{{^3^3)52^3)3^) 


Q5 = 


ivvWv") 


Ql5 


= {{{V3<-P3)4<-P3)3<-P) 


Q6 = 


((^^')3(W")3) 


Ql6 


= f (<^5(W)5) 


Q7 = 


((^V")3(/^")3) 


Ql7 


= f (^5(^V")5) 


Q8 = 




Ql8 


= 


Q9 = 


((^V")5(/^")5) 


Ql9 




Qio = 


= ((W)5(W)5) 


Q20 


= ^'^{^3^3) ■ 



(47) 



the NLO contribution to the flavon superpotential is 

Swi + 6W2 + SW3 



where 



6wi 



12 



Qi 



19 



5W2 



=13 



6w^ 



X20- 



Q 



20 



A 



(4J 



(49) 



It is useful to deal with the contribution Wf + 6wi first. This includes all quartic operators 
that depend on the fields (p, if' and y?" only. In this case the minima of (p, (p' and (p" can 
be analyzed in a analytic form. This part of the superpotential breaks the invariance of 
eq. (HSjl . but is still invariant under the linear trasformation 



ip — )■ QR23 Q ^ip , (p' ^ y fi-R23 ^ 
with Q given in eq. (jHj) and 



p>" -> - QR23 

y 



(50) 



R 



23 




(51) 



with a complex. We start by looking for a minimum for p), pJ and pJ' with the same 
orientation of the one in eq. ( H5l) 



(m,0,0) A 
(0,u',0) A 
(0,0,n")A , 

Along this direction the minimum conditions reduce to 



(52) 



gv!v!' + + 4(x2 + xii)m^ + i(4x6 



4X7 + SXg + 2Xi2)uu'u" = 



M2 1, 
+ — + -(4X6 



4x7 + Sxg + 2xi2)m^ H — (4x4 + 4x8 + 2^10 + Qxi3)u'u" = 
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We find that the values of the u, u' and u" components are the ones given in eq. f HS]) plus 
small perturbations of order VEV/ K. The solution f l52|) is not isolated. It is continuously 
connected to an infinite set of solutions given by 

V9 = fii?23 ^"Vo , = y ^-^23 ^~ Vo ) V^" = -^-^23 ^~Vo • (53) 

y 

Thus the degeneracy present at the LO has been only partially removed by the NLO 
contribution 5wi. The remaining degeneracy is removed by the contribution 5w2- We 
have analyzed the full NLO superpotential Wf + 5wi + 5w2 + Sw-^ by looking for numerical 
solution to the minimum equations. Looking for minima for the fields 99' and if" we 
have frozen the value oi ips to its LO minimum, eq. flTB]) . Under this condition it is easy 
to see that the operators Q13, Qn and Q15 vanish. Moreover the operator Q20 does not 
influence the minima of (p, if' and ip" and the effect of the operators Qis and Qig can be 
absorbed in a redefinition of the parameters Mi and M2 of the LO superpotential. In our 
numerical simulation g, xi-12 and Xi6,i7 are complex random numbers generated with a 
fiat distribution in the square defined by the corners [—(1 + i)/\/2, (1 + i)/\/2]. To get 
an expansion parameter VEV/A of order 0.01, we have taken values of M12/A, {(fs)/-^ 
and A in the square defined by the corners [—(1 + i)/ V2, (1 + i)/ \/2] x 10~^. We have 
performed 50.000 independent minimizations of the scalar potential. We find that the 
mean values of the VEVs, normalized to one up to terms of order (VEV/A)^, are 

(if) = (100, 0.35 + O.m, 0.35 + 0.17i) x 10"^ 

(^') = (-0.51 - 0.03i, 100, -0.53 - 0.05i) x 10"^ 

(y,") = (-0.51 -0.03i, -0.53 -0.05i, 100) X 10^2 _ (54) 

Notice that the induced perturbations are not independent. We have {ip2) = (v^s), (v^'i) = 
{ip'D and (v^g) = {ip2). This is true not only on average, but also separately for each 
individual minimization. The direction of the minima in flavour space is now completely 
determined and coincides, up to corrections of relative order VEV/A with the alignment 
f|T6|) needed to enforce the desired mixing pattern. The only remaining fiat direction is 
that related to the overall scale of and if" (the parameter y in eq. fH3|l ). since also 
the NLO superpotential only depends on the combination <f'(p". This last fiat direction 
is removed at NNLO order where terms depending separately on ip' and cp" first occur in 
the superpotential. Finally the contribution 6w2 also modifies the VEVs of ^, ^' and ips 
compared to their LO values. Also these corrections are of relative order VEV/A. 

In summary the analysis of the scalar potential of the model in the SUSY limit shows 
that the minima of the flavon fields are given by 





[2 VEV'^ 
( Y3*^^ + ^)' P^I^I'P) ^ + '^( ^ 


(0 = 


VEV'^ 
sA + 0{^^) 


{^) = 


VEV^ 

(«,0,0) A + 0(^^) 


{^') = 


VEV'^ 

(0,<0) A + 0(^^) 
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{^") = (0,0,«") A + 0(^^) . (55) 

This proves that the lepton mixing pattern originates from the dynamics of our model and 
not from an ad hoc choice of the underlying parameters. 

6.2 Other higher-order operators 

Beyond the operators (147|) . that correct the lepton mass spectrum through the flavon 
VEVs, there are other higher-dimensional operators contributing directly to lepton masses. 
We only consider NLO contributions. At this order the charged lepton mass matrix mi 
is not affected. At LO mi is dominated by operators of order 1/A. At NLO we find the 
following three invariant operators: 

-l/x^(^'VO^d, 

l^T^{^^'l)H,. (56) 

From the multiplication rules reported in the Appendix and the alignment shown in the 
previous sections it is easy to show that the NLO VEVs of the new triplets are 

(^V')3=(«' u",0,0)A', 

(¥.'V)3 = (0,0,m"m)A2 (57) 

and that they exactly align in the same direction as and respectively. Thus these 
operators do not modify the LO structure of m/. We conclude that, at NLO, the charged 
lepton mass matrix is only modified by the corrections to vacuum alignment analyzed in 
section 6.1. 

The neutrino sector receives corrections from the operators: 

(Wsi^'iy')- (58) 

The first one modifies non-trivially the Dirac neutrino mass matrix mo- The second one, 
after the breaking of the flavour symmetry, can be absorbed by redefining the coupling 
constant yi. The third one changes the Majorana mass matrix M for the heavy neutrinos. 
The neutrino mass matrix m,^ receives two type of corrections at the NLO. One coming 
from the modified vacuum for the fiavon (fs and another one from the operators ( 158|) . 
The corrections to the entries of (ps are unrelated to each other and consequently slightly 
modify the vacuum alignment shown in eq. ( 138|) . Neutrino masses and mixing angles are 
modified by terms of relative order VEV/A. The size of this correction is constrained by 
the agreement between the predicted and observed value of 6*12. Not to spoil the successful 
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prediction of 6*12, the ratio VEV/K should not exceed a few percent. In figure [2] we show 
the relation between sin^ 612 and sin^ 6*13 at the NLO order as a result of a numerical 
simulation with random parameters. The simulation takes into account all the corrections 
coming from eqs. ( l55l) and ( !58l) . 
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Figure 2: Correlation between sin^ 612 and sin^ 6*13 at the NLO. For the triplets ip' and ip" 
the mean values in flMl) were chosen. {ips)/A, and (^')/A have values in the square 

defined by the corners [—(1 + i)10~'^/ \/2, (1 + z)10^^/-\/2]. Similarly the other parameters 
are complex random numbers generated with a flat distribution in the square defined by 
the corners [—(1 + i)/ a/2, (l + i)/ a/2]. The LO prediction for the solar angle is also shown 
(dashed line), along with the extremes of the 3a experimental range (dot-dashed lines). 

We recognize a possible correlation: in general either 612 or 613 can deviate significantly 
from the LO prediction, but not both at the same time. When VEV/A is of order 0.01, 
the maximal correction to the reactor angle remains far below the sensitivity of future 
experiments ^20j and the result from the global fit on neutrino oscillation given in [21], 
where a value not too far from the current limit is reported. On the other hand, 612 could 
be shifted toward the experimental value shown in ([1]). 

7 Lepton mixing angles and Platonic Solids 

On a pure phenomenological basis there are attractive patterns that could provide a good 
LO approximation for the lepton mixing. In particular the data are firmly indicating that 
the atmospheric mixing angle is close to maximal and that the reactor angle is the smallest 
one so that, in a crude approximation, we can take sin 023 = 1/2 and sin^ ^13 = 0. In 
the same approximation, several choices have been suggested for the solar angle, such 
as sin^^ia = 1/2 (BM), sin^ ^12 = 1/3 (TB) and tan ^12 = 1/0 (GR), the choice of the 
present paper. Clearly, depending on the pattern chosen as first approximation to the 
data, appropriate subleading corrections are needed, smaller for the TB and GR mixing 
patterns and larger for the BM one. In the basis where the charged leptons are diagonal the 
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most general neutrino mass matrices corresponding to the considered LO approximations 
are given by 

/x y y\ 

my=\yzw\ , (59) 
\y w z I 

with coefficients x, z and w satisfying the following relations: 

z + w = X BM 
z + w = X + y TB 

z + w = x-V2y GR . (60) 

Thus we have three textures each depending on three independent complex parameters. 
Each of them can be completely characterized by a simple symmetry requirement, that of 
being invariant under two independent commuting parity transformations, U and S: 

S^ = U^ = l ^ [S,U] = . (61) 

Indeed the mass matrix of eq. ( 159|) is the most general one invariant under the so-called 
fi — T exchange symmetry, generated by: 



t/ = 1 . (62) 




Moreover, depending on the particular chosen texture, we have another generator S, which 
can be found by using the constraints in eq. ( !60|l . We list the solutions in Table 4. The 
requirement of invariance under U and S completely determine the mass textures listed 
above. Namely, given a generic neutrino mass matrix m,y the most general solutions to 
the equations: 

U'^ rriy U = rrijy , S'^ rriy S = rriy , (63) 

with U and S given in eq. (l62l) and in Table 4, are the mass matrices defined by eqs. (!59|) 
and (|60|). 

At the same time the requirement that the combination m\mi is diagonal can be fulfilled 
by asking that the charged lepton sector is invariant under the action of a cyclic symmetry 
of Zn type. Calling T the corresponding generator, with T" = 1, we can choose T so that 
the solution of 

{m\mi) T = (m|mi) (64) 

is a diagonal matrix. In minimal constructions realizing BM, TB and GR mixings the 
generator T can be chosen as in Table 4. With this choice T has the additional property: 

{STf = 1 . (65) 



A model giving rise to the desired lepton mixing matrix Upmns can be obtained starting 
from the family group Gf generated by S, T and U. The family symmetry should be 
spontaneously broken by a set of scalar fields in such a way that, at the LO, the charged 
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Table 4: Generators S and T for the different mixing patterns (0 = (1 + V5)/2 and 
CO = -l/2 + iV3/2). 



lepton sector has a residual invariance under the group generated by T, whereas the 
neutrino sector has a residual invariance under the group generated by S and U. The 
desired lepton mixing Upnms arises by construction, independently from the base choice. 
In existing models the — r exchange symmetry generated by U arises at the LO as 
an accidental symmetry. In this case the family symmetry G/ is generated by S and T 
satisfying: 

= {STf = 1 and T" = 1 , (66) 

with n = 4 for BM, n = 3 for TB and n = 5 for GR. These are the defining relations 
of 5*4 (BM) [5], (TB) [22] and (GR), which are the proper symmetry group of 
the cube/octahedron, tetrahedron and dodecahedron/icosahedron, respectively. We find 
rather intriguing that Platonic solids and their symmetries are in a natural correspondence 
with the most popular lepton mixing patterns. Notice that these groups are all subgroups 
of the modular group, defined by the presentation S'^ = {ST)^ = 1. 

8 Conclusion 

We think that the GR mixing pattern, where sin^^23 = 1/2, O13 = and tan ^12 = 1/0, 
should be considered on the same foot as other more popular schemes, such as the TB and 
the BM ones, in our attempts to construct a model of lepton masses and mixing angles. 
Indeed the GR scheme is compatible with the experimental data. The largest deviation is 
for the solar angle, where the predicted value is about two standard deviations below the 
present experimental central value. In this work we have built a supersymmetric model 
reproducing the GR pattern in the LO approximation. In the limit of exact GR mixing we 
have identified two transformations S and T, leaving invariant the neutrino mass matrix 
and the charged lepton mass matrix respectively, and generating the discrete group A^. 
Following this hint, we have chosen as the family symmetry of our model A^ x Z5 x Z3 
where the Z5 x Z3 factor forbids unwanted couplings between the fiavon fields and the 
matter fields. In the supersymmetric limit we have analyzed the most general scalar 
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potential for the flavon fields up to terms suppressed by one power of the cutoff A. In a 
finite portion of the parameter space, without any fine-tuning of the parameters, we find 
an isolated minimum of the scalar potential where the fiavon VEVs give rise to the GR 
mixing pattern, up to terms of order VEV/K. Choosing VEV/h. of order few percent 
we can have an excellent agreement between theory and data for both the solar and the 
atmospheric mixing angles. The mixing angle ^13 is expected to be of order few degrees. 
The neutrino masses depend on three complex parameters so that both types of ordering 
can be accommodated. Neutrino masses and squared mass differences can be fitted but 
not predicted. At the LO we find restrictions on the allowed value of rngg, once the mass 
ordering and the smallest neutrino mass have been fixed. 

To achieve the desired vacuum alignment there is no need of driving fields, a tool often 
used in this type of constructions. Neglecting matter multiplets, the energy density of the 
theory only depends on the fiavon fields, that are self-aligned at the minimum. To our 
best knowledge our model is the first example where the GR mixing pattern is derived 
from a full minimization of the energy density. 

We think that our model provides a valuable alternative to other existing proposals. 
There are not sufficient hints in the data to prefer other mixing patterns, such as the TB 
one, to the GR one. Indeed in most of the existing models, including the present one, 
TB, BM, GR or other mixing patterns are only lowest order approximations, unavoidably 
corrected by powers of the symmetry breaking parameters. From this point of view, 
TB and GR mixing patterns can be both considered excellent first order approximations 
to the existing data. It is remarkable that the TB, BM and GR mixing patterns can 
be obtained from minimal constructions based on the symmetry groups A4, 5*4 and A5, 
which are the proper symmetry groups of the Platonic solids. To make a comparative 
experimental test of these constructions other observable quantities should be considered, 
such as for instance the rates of lepton flavour violating processes, which, depending on 
the assumed supersymmetry breaking scale, could be within the reach of the presently 
running or planned experiments. 
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Appendix A. Kronecker products 

We report here the complete hst of the Kronecker products for the group A^. We as- 
signe a — (ai, 02, 03)'^ and h — (61, 62, ^'3)'^ to the 3 representation, while a' — (a'^^, Og, a'^^ 
and h' — {h'nh^ih'^^ belong to the 3' representation, c = (ci, C2, C3, C4, 05)-^ and d — 
{di, d2, d3, dA, d^Y are pentaplets; / = (/i, /2, fif and g = (t/i, g'2, t/a, QaY are tetraplets. 
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3 ® 3 = 3„ + (1 + 5), 



1 = aibi + 0263 + a3^'2 

3 = (0263 - 03^*2, ai&2 - 02^1, oa&i - oi&a)^ 

5 = (a^^^ _ _ ^(a^52 + a2&i), -^^02^2, -y^a3&3, --^(ai&3 + 0361))^ 

3' ® 3' = 3'a + (1 + 5), 

1 = a[b[ + 0/263 + 0,362 

3' = (0263 - 462, a[b2 - a'^b'^, a'^b\ - a[b'^f 



5 = {a[b[ 



3' = 4 + 5 



2 ' 



"'^^^ "^^'^ ^^a',b',,-y^{a[b', + a',b[),-^ia[b', + a^^ 



A ( iJ ^3^*2 ,/ , 0^363 ,/ 02^*2 u , <^2b3-.T 

_ , 0261 + 1/20362 O162 + V2O363 O163 + A/2O262 O3 61 + 1/20263 J, 

' ^ 7! — ' 71 — ' 7! — ' 7! — 

3(g)4 =3' + 4 + 5 

1 1 
3' = (a2g4- a3gl,--^{^/2alg2 + a2gl + asg3),--^{^/2algs + a2g2 + asg4))'^ 

4 = {aigi + i/2o35f2, -aifi'2 + 72o25'i, ai5'3 - 72a35'4, -ai5'4 - \^a2g^f 

5 = (035^1 + 025^4, (Veals'! - a35'2) , (72ai5'2 - 3035^3 + a25'i) , 

-^{V^aig^ - 3025-2 + a^g^), ^j^i-^aig^ + 025-3))^ 
3' 04 = 3 + 4 + 5 

1 1 

3 = (Gas's - 4fl'2, -^^ol^9\ + «2fl'4 - ass's) ! (-72a'ig'4 + a2fl'2 - oUi))^ 

4 = (o'lg-i + a/2o35(3, o'i5f2 - 72a3fi'4, -a'lSs + V^Oafi-i, -a'iS4 - \/2o2fi'2)^ 

5 = (03^2 + 02^3, ^{^Pid^gx - 30254 - 4^3) , - y|(72oi52 + 4^4) , 

-y|(72o'i53 + 4^1), ;^(-72o'i54 + 30351 + 0252))^ 
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5 = 3 + 3' + 4 + 5 

3 = (^^^ + a3C2 - a2C5, + aiC2 - v^oacs, - aiC5 - "\/2a2C4)^ 
. . , a2C5 - a3C2 aiC3 + v^(a3C4 - 0202) aiC4 + -y2(a2C3 + 0305) ^ 

4 = (4aiC2 + 2"\/3a2Ci + "\/2a3C3, 2aiC3 - 2\/2a2C2 — SV^a^Ci, 

2aiC4 — 3-\/2a2C3 + 2-\/2a3C5, —AaiC^ + -\/2a2C4 + 2-\/3a3Ci)^ 
aiC2 + "\/2a3C3 2aiC3 + V^a2C2 



(0205 + 0302, a2Ci - 



3' 05 = 3 + 3' + 4 + 5 

, 0303 + a2C4 -a[c2 + \/2(a3C4 + ^cs) a'^cs + \/2(a2C3 - 0302) , ^ 

3 = („,e, + _^,_ ,_ ) 

3' = ~ "3'^^ ~ "2'^^' ~ ^'^^^ ~ V2a'2C5, - «iC4 + \/2a2C2)^ 

4 = (2aiC2 + 3^2405 - 2V2a^C4, -Aa'^Cs + 2^3401 + v^agCg, 

— 4aic4 — V2a2C2 + 2-\/3ci3Ci, — 2a'^C5 — 1\pld.f.-i, — ?>\pi(i^C2)^ 



5 = (a2C4 - 0303, , -a2Ci , 

/ a'^C4 + -\/2a2C2 — 20^05 + -\/2a2C3xT 

+ — 7i — ' ^ 

4 ® 4 = (3 + 3')a + (1 + 4 + 5), 

1 = /l5'4 + /2fi'3 + + /45'1 

3 = {fi94- f49i + f 392 - f 293, V2{f 294 - 1492), V2{fig3- fsQi))^ 

3' = ifi94- f49i + f293 - f392,V2{f3g4- f4g3),V2{fig2- f29i)V 

4 = (/aS's - /45'2 - f294, fi9i + /35'4 + f493, -f494 - fi92 - f29i, -f292 + figs + f39l) 



T 



ifm + f4gi - hg2 - /2^3, -y ^(2/3^3 + /2«/4 + 7452), ]J'^{-2figi + 73^4 + /4^3), 
^(-2/4^4 + /2^i + /i^2), y|(2/2^2 + figs + f3gi)f 



4(8)5 = 3 + 3' + 4 + 5 + 5 

3 = (4/iC5-4/4C2-2/3C3-2/2C4,-2v^/lCi- V2(2/2C5-3/3C4 + /4C3), 

V2{-fiCi + 3/2C3 + 2/3C2) - 2v^/4Ci)^ 

3' = (2/iC5-2/4C2 + 4/3C3 + 4/2C4,-2v^/2Ci + ^(2/4C4 + 3/iC2-/3C5), 
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V2(/2C2 - 3/4C5 + 2/1 C3) - 2v^/3Ci)^ 
4 = (3/lCi + V^(/2C5 + /3C4-2/4C3),-3/2Ci + y6(/4C4-/iC2 + 2/3C5), 
-3/3C1 + a/6(/iC3 + /4C5 - 2/2C2), 3/4C1 + a/6(/2C3 - /3C2 - 2/iC4))'^ 

51 = (/1C5 + 2/2C4 - 2/3C3 + /4C2, -2/iCi + V6/2C5, /2C1 + y^(-/iC2 - /3C5 + 2/4C4), 

-/3C1 - (/2C2 + /4C5 + 2/1C3), -2/4C1 - V6/3C2)^ 

52 = (/2C4 - /3C3, -/iCi + ^^^^^ ^,-y|(/lC2 + /3C5-/4C4), 

/2/. , . , / N / 2/3C2 + /lC4 + /2C3^ 

-y gl/ica + /2C2 + /4C5), -/4C1 ) 

5 (8) 5 = (3 + 3' + 4)„ + (1 + 4 + 5 + 5), 

3 = (2(04(^3 - 03(^4) + C2(i5 - C5(i2, V3(c2«il - Ci(i2) + \/2(c3(i5 - 05^3), 

Vsic^di - cid^) + v/2(c4(i2 - 02^4))'^ 

3' = (2(c2(i5 - 05(^2) + C3(i4 - 04(^3, V3(c3(ii - Ci(i3) + \/2(c4(i5 - 05^4), 

A/3(ci(i4 - C4(ii) + A/2(c3(i2 - £2^3))^ 

(C3(i5 + 054) - 4c4(i4 / , , , X (C44 + £5(^4) - 4c2(i2 

4s = ((Ci(i2 + C20(i) ,-(Cl(i3 + C30(i) , 

/ J , J\ (02^3 + C34) + 4c5(i5 . , , (C2d4 + 044) + 403^3 ;r 
(CiO(4 + C4di) -j^ , (Ci4 + Cgdi) --^= ) 

4a = ((Ci4 - C24) + y/^(c34 - C54), (Cl4 - C34) + ^^(044 - C54), 
(04^1 - Ci4) + y|(c34 - C24), (Cl4 - C5(il) + y^(c44 - C24))^ 

5i = (ci4 + C24 + C54 + ^^^^ ^ ^^^^ , -(ci4 + C2di) + y^C44, ^(ci4 + C34 - v^(c44 + C54)), 

^(ci4 + C44 + "\/6(c24 + C34)), -(ci4 + C54) - 'y|c34)^ 

_ 2ci(ii + C24 + C54 -3(ci4 + C2(ii) + \/6(2c44 + C34 + C54) 2c44 + 2c54 + C24 
' " ^ 2 ' 6 ' 7! ' 

2c24 + 2c34 - C54 -3(ci4 + c^di) + \/6(-2c34 + C24 + C44) xt 

1 = Cidi + C34 + C44 — C2rf5 — C54 
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Appendix B. Other representations of A 



Shirai's basis. In this section we show how the basis chosen in this work is related to 
the representations given in previous articles on the group A^. In [HI |23] the Shirai base 
was used, in which the presentation of the group is given in terms of two matrices Ssh and 
Tsh satisfying the following algebra: 

^sh = = (Tgf^SshTshSshTgi^SshTshSshTgil)^ = 1 (67) 

In order to connect them to the matrix S and T satisfying 

= T^ = {STf = 1 (68) 

with T diagonal, we take an intermediate step and define first 

S' = Ssh (69) 

and a matrix T' such that {S'T'Y = 1 Then, we define Tgf^SshTgf^SshTg,^ SshTshSshT^h = 
Ash- Since S'T' = Ash we obtain 



r = s'-^A 



Sh 



S'A 



Sh 



SshAsh 



(70) 



Note that T' is not diagonal. To make it diagonal, we let an unitary matrix U act on it, 
such that 



and finally 



S = U^S'U. T = U^T'U. 



S = U^SshU, T = U^S.,hA,^,U. 



Sh^Sh^ 



(71) 



(72) 



(73) 



Cummins-Patera's basis. The Cummins-Patera's basis [HI [19] is generated by two 
elements with presentation 

^2 = 53 = (AS)5 = 1, 

and the quintic representation are explicitly given by 
/I ON^ 

1 

~2 



A 



-1 





Vo 




-10 
1 
001/ 



B 



/ 


1 





1 




-UJ\ 







1 


1 


-1 


-1 




-1 


-1 











u 


1 







—u 


\ 




1 






/ 



(74) 



where u = e"^ is the cubic root of unity. We can define a unitary transformation Ucp 
that relates the elements of the Cummins-Patera basis to the one introduced in Section 3 
as follows: 

S = UlpAUcP, (ST) = UlpBUcP, T = Ulp{AB)Ucp. (75) 
The unitary matrix is 

/ 

-\/3a;3 



CP 



1 

7E 



_i 

U 2 



-u 



tu 2 
2 



1 

2 

ia 



U) 2 , 



— e 



-17 



2 

1 

20 



_l 
-00 2 



2 



UJ 2 . 



24, ' 



_ 1 

UJ 2 , 



1 

24, 
ia 



—UJ 2 



— e 



1 

2 



(76) 



24 



where a = arctan (— /3 = arctan ( y 3 — and 7 = arctan ( y |). It is straightfor- 
ward to verify that acting with U^p on the vector (0, 0, 0, Z, Z)'^ gives a new vector with 
the form given in eq. (138|) . 
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